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ABSTRACT 


In  this  paper  the  problem  of  fault  analysis  in  Linear  and  Affine 
Sequential  Circuits  is  treated.  These  classes  of  systems  provide  for 
the  treatment  of  several  linear  and  nonlinear  faults  common  in  digital 
circuitry.  The  solution  to  the  linear  and  affine  sequential  circuit 
fault  analysis  problem  is  obtained  via  the  development  of  a spectral 
theory  for  such  systems  over  finite  fields.  A stepwise  fault  analysis 


procedure  for  this  problem  class  is  presented  along  with  many  examples 


illustrating  the  advantages  that  memory  provides  in  digital  fault  analy- 
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CHAPTER  1 


I 


INTRODUCTION 

1.1  Introduction  to  Fault  Analysis: 

The  rapid  progress  of  modern  integrated  circuit  technology  has 
made  it  possible  to  manufacture  many  new  digital  circuits  with  complex- 
ity orders  of  magnitude  greater  than  the  circuits  of  the  preintegration 
period.  A semiconductor  "chip"  which  is  just  the  size  of  a fingernail, 
may  have  many  individual  circuit  elements  embeded  in  it.  Moreover,  a 
number  of  such  chips  may  be  combined  together  to  construct  a large  cir- 
cuit. This  renders  the  structure  of  the  circuit  more  complex  which  in 
turn  makes  maintenance  difficult.  If  the  circuit  does  not  behave  in  the 
manner  it  is  supposed  to,  the  faulty  plement  in  the  circuit  must  be 
located.  Techniques  used  for  locating  the  faulty  part  or  parts  of  a 
circuit  are  termed  as  fault  analysis  techniques. 

Historically,  fault  analysis  techniques  for  analog  and  digital 

circuits  have  developed  independently.  The  reader  is  referred  to  the 

I 2 

bibliographies  compiled  by  Rault  ’ for  a review  of  the  literature  in 
the  two  areas.  Digital  fault  analysis  techniques  are  mainly  combina- 
torial in  nature,  i.e.,  ore  *‘'rts  each  component  of  the  circuit  by 
applying  a family  of  test  inputs  to  that  circuit.  If  one  obtains  the 
expected  response  from  the  circuit  for  each  test  input,  the  circuit  is 
operating  properly  and  so  are  all  of  its  components.  On  the  other 
hand,  if  the  circuit  fails  to  operate  correctly  for  one  or  more  test  in- 
puts, the  faulty  component(s)  may  he  isolated  by  determining  the  s- * ..f 
components  which  are  exercised  by  exactly  that  set  of  input  test  sig- 
nals. In  the  case  of  analog  circuits,  one  has  the  advantage  of  testing 

the  circuit  by  exponential  signals  of  different  frequencies.  Hence  the 

1 
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gain  between  various  circuit  test  points  at  each  of  several  different 
frequencies  may  be  measured.  Then  by  some  means,  an  equation  solver  or 
optimization  routine,  the  circuit  parameters  are  determined  which  will 
yield  these  gains.  If  the  resultant  component  parameters  are  in  the 
operating  range,  the  corresponding  component  can  be  assumed  to  be  work- 
ing correctly,  whereas  if  a computed  component  parameter  deviates  sig- 
nificantly from  its  operating  range,  the  corresponding  components  may  be 
assumed  faulty.  * 

Fault  analysis  techniques  for  digital  and  analog  circuits  with 
memory  differs  from  each  other.  In  digital  circuits,"*  due  to  combina- 
torial approach,  test  comp'exity  increases  exponentially  with  memory.  J 

(Since  the  number  of  entries  in  a circuit  truth  table  increases  expo- 
nentially.) In  cast  of  dynamical  analog  circuits,  however,  just  the 
opposite  is  true.  One,  here,  uses  tne  mui tit'requency  testing  techniques  > 

in  which  one  may  determine  the  gain  between  a pair  of  test  points  at 
several  different  frequencies  simultaneously  from  a single  test  signal. 

Thus  total  number  of  test  inputs  required  for  testing  dynamical  analog  * 

circuits  are  less  than  for  a meinoryless  circuit  of  a similar  complexity. 

1.2  Fault  Analysis  of  An..l'«a  Circui  ts : ^ ,P' 

Analog  fault  analysis  is  heavily  pedicated  on  spectral  theoretic 

technigues  formulated  in  a frequency  domain  setting.  In  the  area  of 
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fault  analysis  of  analog  circuits  a number  of  recent  papers  ’ ’ have 
given  considerable  attention  to  the  component  connection  model  of  a 
large  scale  system.  Via  the  component  connection  model  fault  anal 
is  based  on  a system  modeling  technique  wherein  the  observable  system 
behavior  is  expressed  explicitly  as  a function  of  the  internal  re- 
sponses. In  the  component  connection  model,  one  assumes  that  the  ith 
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component  is  characterised  by  an  operator  equation  (usually  expressed  in 
the  frequency  domain). 


bi  = Z1  ai  1 - 1 ,2 ,-,-,-,n  (1.1) 

mapping  the  component  input  vector,  a-,  into  its  output  vector,  . 
Although  in  actual  practice  one  normally  works  with  the  n separate 
component  equations,  notationally  (1.1)  may  be  combined  into  the  single 
equation 


b = Z a 


0.2) 


where  b = col . (Ik  ) 
a = col . (a.j ) 

Z = diag  (Z^ 


Since  the  connection  elements  (adders,  scalers,  etc.)  are  all 
algebraic,  the  connection  model  is  represented  by  a set  of  linear 
algebraic  equations  of  the  form. 
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(1.3) 


where  y is  the  overall  system  output  and  u is  the  input  vector  for  the 
overall  system  which  are  related  by  the  overall  system  operator,  S,  by 
y - S u (1.4) 

Upon  combining  equations  (1.2),  (1.3),  and  (1.4),  one  obtains  the 
equation. 


s = i22  * l21(i-zl11)"1zl12  = f(Z)  (l.S) 

relating  the  overall  system  operator  to  the  composite  component  operator. 
Since  the  function,  f,  is  entirely  determined  by  the  fixed  connectior 
matrices,  it  has  been  termed  the  connection  function.  Now,  if  Z and 
S are  characterized  in  the  frequency  domain,  one  obtains. 


S(.o)  = f(Z(«)) 


(1.6) 


4 

since  the  connections  are  independent  of  ^rtquency. 

Finally,  in  many  practical  systems  one  can  write  the  component 
matrix,  Z(l.)  , in  the  form 

Z(u)  = E(u)CFU)  (1.7} 

where  E(ui)  and  F(oj)  are  frequency  fjpendent  matrices  determined  by  the 
component  types  but  not  their  actual  values  and  C is  a frequency  inde- 
pendent matrix  of  component  values  (i.e.  for  an  inductor  the  j.„  goes  in 
E or  F and  L goes  in  the  C matrix).  Such  a viewpoint  is  quite  reasonable 
for  the  fault  analysis  problem  wherein  one  can  reasonably  assume  that  the 
component  types  remain  fixed  and  all  faults  manifest  themselves  as 
changes  in  component  values.  We  then  combine  the  known  E and  F matrices 
into  the  connection  function  and  characterize  our  system  by 

SU)  --  f(E(,,)CF(c  ))  » f (C)  (1.8) 

Equation  (1.8)  is  just  the  right  form  in  which  to  study  tne  fault 
analysis  problem  for  if  we  make  the  (standard)  hypothesis  that  all  faults 
take  the  form  of  errors  in  C with  the  connections  and  component  types 
fixed,  i.e.  f is  fixed,  then  one  merely  measures  S(u)  at  some  freciuencv 

Ci» 

and  solves  equation  (1.8)  for  C.  Unfortunately,  solution  of  these  non- 
linear equations  is  equiv,:l"’t  to  left  invertibi  1 ity  of  the  matrix 

K = l{2  ® l?1  (1.9) 

which  requires  that  the  system  have  a large  number  of  test  points. 

If  equation  (1.8)  is  not  soluble  one  has  two  alternatives;  either 
add  more  test  points  (which  increases  hardware  costs'  or  use  several 
test  frequencies  (which  increases  software  costs).  The  former  is 
straightforward  and  may  be  carried  out  algorithmically  by  adding  addi- 
tional test  points  in  such  a manner  that  the  additional  rows  added  to  the 
matrix  K will  render  its  columns  linearly  independent.  From  a practical 


5 

point  of  view,  however,  it  is  preferable  to  use  additional  test  frequen- 


cies at  the  same  test  points  in  which  case  one  must  solve  the  set  of 
simultaneous  equations 

SUV,  = f U) 


su) » f ft: 

C u-'*' 

c 


' ' V " t(,1  -r  • 


(1.10) 
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ity.  Since  a linear  Sequential  Circuit  can  be  viewed  as  a linear 
operator  on  a sequence  space,  a perfectly  valid  spectral  theory  for  the 
Linear  Sequential  Circuit  may  be  formulated.  The  resulting  spectral 
theory  parallels  the  steady  state  frequency  domain  theory  for  analog 
circuits  and  thus  may  be  used  to  formulate  a fault  analysis  procedure  for 
Linear  Sequential  Circuits  which  closely  parallel  the  multi  frequency 
testing  techniques  for  analog  circuits. 

The  sequential  circuits  are  defined  over  a finite  field.  These 
finite  nelds  are  denoted  by  GF(p)  where  p is  0 prime  and  its  nth  order 
extension  is  denoted  by  GF(pn).  A detailed  review  of  finite  fields  and 

an  algorithm  for  generating  elements  of  GF'ofl)  is  given  in  Appendix  A. 

13 

All  these  extension  fields  lie  in  the  algebraic  closure  of  a finite 
field. 

in  the  Chapter  H,  fr.e  pectral  tneory  for  Linear  Sequential  Cir- 
cuits is  formulated.  For  this  purpose  a Linear  Sequential  Circuit  is 
mathematical ly  described  by  a pair  of  difference  equations  over  a finite 
field.  Rather  than  interpreting  this  set  of  equatiors  as  the  tradition- 
al initial  value  problem,  it  is  interpreted  as  a central  value  problem 
of  finding  two  sided  state  output  sequences,  i.e.,  for  the  posi- 
tive and  negative  values  of  timp. 

The  delays  used  in  the  sequentia1  circuit  are  interpreted  as 
predictors,  i.e..  The  next  state  value  of  the  sequential  circuit  is 
the  delayed  present  stats,  of  the  sequential  circuit.  The  corresponding 
definition  of  the  D-t-^nsform  is  given  in  Append i x B. 

The  entire  mathematics  outlined  in  the  previous  section  for  the 
analog  case  goes  through  if  one  interprets  the  vectors;  a,  b,  u and  y; 
as  sequences  taking  their  values  in  the  finice  field  rather  than  real 
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valued  functions  and  assumes  that  the  entries  in  the  L matrices  tatce 
their  value  in  the  same  finite  field. 

In  the  case  of  sequential  circuits  if  one  deals  with  D-transform 
rather  than  laplace  transform,  equation  (1.1)  through  equation  (1.6)  goes 
through  i .e. , 

S(D)  = f (Z(D ) ) = f ( (G)F(D) ) = fD(G)  (1.11) 

where  f^  is  a nonlinear  function  which  is  entirely  determined  by  com- 
ponent dynamics  and  the  connections  relating  component  parameter  values 
to  the  systems  input  and  output.  Here,  the  implication  that  is 
known  and  non-faulty  is  that  all  faul ls  occur  in  the  scalers,  G,  with 
memory  elements  and  connection  good.  In  particular,  this  implies  a 
l.'eac  system  fails  linearly  and  hence  one  may  include  "stuck-on-zero" , 
"open"  and  "short  circuit"  faults  but  not  "stuck-on-one"  faults  which 
are  nonlinear.  "Stuck-on-one"  faults  are,  however,  included  when  a 
Linear  Sequential  Ci?'cuit  is  generalized  to  the  case  of  affine  circuits 
which  fail  affinely.  The  latter  generalization  also  permits  a slight 
general ization  of  the  traditional  Linear  Sequential  Circuit  case  by 
allowing  NOT  gates  and  bias  sources  in  addition  to  the  usual  Linear 
Sequential  Circuit  components. 

Finally,  the  analog  of  exponential  test  functions  for  the  case  of 
sequential  circuits  are  the  sequences,  {ue(,  of  the  form, 

luke}  = { e k } , k = 0,41  

where  c is  an  element  of  the  algebraic  closure  of  a finite  field. 

The  input  sequence  { u ^ } = {e  } yields  the  equality, 

yk  = S(e)  ek 

for  an  appropriate  initial  state. 

These  exponential  sequences  "live"  in  an  extension  of  the  space  on 
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which  the  actual  circuit  is  defined  and  therefore  can  not  be  physically 
implemented.  This  is  precisely  the  same  phenomena  which  occurs  in 
analog  circuits  where  one  must  test  with  the  real  valued  input  Sin(u:t) 
to  obtain  information  about  the  complex  valued  eJut.  Therefore,  in 
Chapter  II,  a technique  for  computing  S(e)  di-ectly  from  the  system 
impulse  response  is  derived. 

Once  S(e)  is  known,  one  can  write  the  set  of  equations  parallel  to 
those  used  for  analog  fault  analysis  for  several  e's  in  the  algebraic 
closure  of  a finite  field,  i.e.. 


S{e}) 


fe  (G) 
el 


S(e2)  - fe?(G) 


(1.12) 


As  in  the  analog  case  it  may  be  possible  to  solve  this  set  of 
simultaneous  equations  even  though  no  single  equation  has  a unique 
solution.  Thus  it  is  possible  to  exploit  the  dynamics  in  a Linear  Se- 
quential Circuit  in  a similar  manner  to  that  used  in  the  analog  case,  so 
as  to  simplify  the  fault  analysis  procedure. 

In  Chapter  II,  the  required  spectral  theory  is  formulated  and  also 
a technique  for  computing  S ( e ) directly  from  the  impulse  response  is 
derived.  In  Chapter  III,  a formula  for  computing  S(e)  from  the  linear 
Sequential  Circuit  component  parameters  and  the  connection  matrices  is 
obtained.  As  such,  one  may  test  a circuit  with  an  impulsive  input, 
compute  S(e)  therefrom  and  then  compute  the  component  parameter,  by 
inverting  this  latter  formula.  This  inversion  process  can  be  formulated 
as  the  solution  of  a set  of  polynomial  equations  in  several  variables 


via  the  "term  expansion  algorithm"^  which  is  described  in  Chapter  III. 

A family  of  illustrative  examples  appears  in  Chapter  IV  and  the  general- 
ization to  affine  circuits  is  presented  in  Chapter  V. 

* 

I 

ft 

l 

ft 

ft 

ft 
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CHAPTER  II 


SPECTRAL  THEORY  - LINEAR  SEQUENTIAL  CIRCUITS 
2.1  "Two  sided"  Linear  Sequential  Circuits: 

1 5 

The  characteristics  of  a Linear  Sequential  Circuit  (LSC)  are 
depicted  in  Figurt  21. 


1 


Input 


Figure  2.1:  Block  Diagram  Representation  of  a Linear  Sequential  Circuit 


Essentially,  a Linear  Sequential  Circuit  consists  of  storage  ele- 
ments and  combinatorial  logic.  An  input  sequence  applied  to  an  LSC 
results  in  an  output  sequence,  whose  present  value  is  a linear  function 
of  the  present  input  value  and  the  present  state.  The  present  state  in 
turn  is  a linear  function  of  past  states  and  past  inputs.  Hence,  a 
Linear  Sequential  Circuit  can  be  viewed  as  a linear  operator  on  a 
sequence  space. 

Mathematical ly , a "two  sided"  Linear  Sequential  Circuit  over  a 
finite  field  is  represented  by  a set  of  difference  equations. 
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xk+1 


Axk  + 3uk 
Cxk  + Juk 


- ® < k < - (2.1 


Where  the  state  sequence  X,  the  output  sequence  Y and  the  input 
sequence  U are  respectively  s,  m and  r dimensional  vectors  over  the 
finite  field.  xfc,  yk  and  uk  are  respectively  the  state,  the  output  and 
the  input  sequence  values  at  time  k for  the  LSC.  The  matrices  A,  B,  C 
and  J are  finite  field  valued  and  constant  with  dimensions  consistent 
with  xk,  yk  and  uk- 

Traditionally  one  interprets  equation  (2.1)  as  an  initial  value 
problem  in  which  one  desires  to  find  "one-sided"  sequences;  xk  and  yk, 
k 0 satisfying  equation  (2.1)  for  a given  sequence  uk  and  initial 
state,  = x^.  One  in  arpretes  equation  (2.1)  as  a central  value  prob- 
lem wherein  one  seeks  "two  sided"  sequences;  xk  and  yk,  - • < k < ® 
satisfying  (C.l)  for  a given  sequence  uk  and  central  value  xQ  = x. 
Unlike  the  case  for  LSC's  defined  over  "one  sided"  sequences,  this 
central  value  problem  for  difference  equations  (2.1)  does  not  admit  a 

unique  solution  for  all  x.  = x.  In  this  chapter,  the  theory  of  exist- 

0 

ence  and  uniqueness  of  solutions  to  the  central  value  problem  is  de- 
veloped and  a viable  spectral  theory  for  the  difference  equations  (2.1) 
which  closely  parallels  steady  state  frequency  domain  theory  for  con- 
tinuous linear  systems  is  formulated.  Such  a spectral  theory  provides 
a way  for  determining  faults  in  a Linear  Sequential  Circuit. 

The  existence  and  uniqueness  of  solutions  to  the  homogeneous 


version  of  (2.1)  which  is. 
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xk+l  = Axk 

yk  = Cxk  - - < k < co  (2.2) 

xo  = S- 

will  now  be  presented  in  the  first  two  lemmas.  The  existence  and 
uniqueness  of  solutions  to  (2.1)  w’ll  then  be  established  in  the  third 
1 ermia . 

Lenina  1:  Let  A be  a linear  transformation  on  a vector  space  X, 

then  X has  a "fitting  decomposition"^  given  by, 

* ■ f„  © f, 

Where  fQ  = {z|Atz  = 0,  z e X,  t > T for  some  T}  and  f1  = R(At)  = 

R(At+1)  « t > T. 

The  subspaces  fQ  ard  f^  are  respectively  the  "Fitting  Null"  and 

the  "Fitting  1"  components  of  X,  Furthermore,  the  restriction  of  A to 

0 

V A0  : fQ  * fQ  is  nil  potent  i.e.,  Aq  =0  for  some  i,  and  the 

restriction  of  A to  : f-j  ► f,  is  invertible  on  f ^ . 

The  proof  of  this  lemma  is  given  in  the  literature,^  but  for  com- 
pleteness and  convenience  is  included  here. 

Proof: 

2 

Let  x f R(A  ),  then  x * A[Ay'J  for  some  y.  Which  implies  that 
x t R(A)  and  hence  R(A^)  £R(A).  It  then  follows  that  rank  (A^)  <_ 

p 

rank  (A).  Similarly,  one  can  inductively  show  that,  R(A)  D R(A  ) s 

R(A^)  d D R(Ar)  and  rank  (Ar)  <_  rank  (Ar  ^ ) <_ < 

rank  (A“)  <_  rank  (A).  Mote  that  as  r gets  larger,  rank  (A  ) will  r vir- 
tually remain  constant  since  it  is  bounded  from  below  by  zero.  Thus 
there  exists  an  r such  that 


i 


I , 


1 * 

: 

? l 
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rank  (Ar)  = rank  (Ar+1)  = or 

R(Ar)  = R(Ar+1)  = » fr 

Define  = (z | A1 z = 0,  z c XI  (2.3A) 

2 

Then  B-j  = { z | Az  = 0 = A z = A[B|]},  so  one  has  B-j  c_  similarly  one 

can  inductively  show  that  Bj  c.  B2  c b3  c Then  there  exists 

s for  which  the  null  space  defined  by  equation  (2.3A)  will  attain  some 

constant  value.  It  follows  that  B = £ = . . . = f 

s s+i  * o’ 

Let  t = max(r,s).  Then  f = Bt  and  f1  = R(At) . 

Any  x e X can  be  written  as  x = (x  - A*y)  + A*y,  Afcy  e f-|  for  some 

y e X.  We  also  know  that  A^x  = A2ty  i.e.,  At(x  - A*y)  = 0 which 

implies  (x  - At y ) t f_.  Hence  we  have  x = f + f , . 

0 0 t 

Let  z c fQ  fl  f-j.  Then  z = A^w  for  w c X,  since  aS*  i f , A*w  = 0 

implying  z = 0 = fQ  fl  f ^ , thus  X = fQ  ® f 1 . 

Since  f^  = R(Ar)  = R(Ar  ^ ) = = r^A,  hence  A is  surjective 

in  f.|.  Since  fQ  = B+,  A*  = 0 in  fQ  then  A is  an  isomorphism  on  f^  i.e. 
A^  : f,  •+  f^  is  invertible  on  f ^ . 

Lemma  2:  (Homogeneous  Case) 

(a)  The  Equation  xk+1  = Axk,  - ®<k<«>,  X()  = x (2.3) 

defined  over  a finite  field  has  a solution  if  and  only  if  x = x E f, . 

o — 1 

In  this  case  the  solution  is  unique,  takes  its  values  in  f^  for  all  k 
and  is  given  by  xk  = Aj  x,  - <*>  < k < (2.4, 

(b)  The  Equation  xfe+1  = Axk>  yk  = Cxk>  - « < k < - (2.5) 

xQ  - x defined  over  a finite  field  has  a solution  if  and  only  if  x i f 1 . 
in  this  case  the  solution  is  unique,  takes  its  values  in  C(f.j)  for  all  k 

|r 

and  is  given  by  yk  = CA-j  x^  - «*>  < k < <*■.  (2.6) 


< k < 


Proo^: 
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(a)  The  equation  x^,  * Ax^,  - «°  < k < ® can  also  be  written  in 


the  form,  = ATxkT  t R(AT)  c 


(2.7) 


The  condition  that  2L  £ ^-|  follows  from  the  fact  that  in  any  solution 


x « xQ  can  be  written  in  the  form, 


x = xQ  = ATx_t  e R(AT)  = f1 

k 

The  claim  is  = A^  x^  xQ  = x.  t f^  (recall  that  the  inverse  of 
exists),  - ” < k < °°  is  a solution  to  (2.3).  lo  prove  this  show  that 
(2.4)  satisfies  (2.3). 

Ax^  - A A^x^ 

= A^  x [A  can  be  replaced  by  A1 

= A]*'+1j(.  since  A^kx  e f^] 

= *k+l 

Uniqueness  of  the  solution  can  be  proven  by  contradiction.  Let 
there  be  two  solutions  e f ^ and  x^  e f ^ . Such  that  x^  = x^  - x.- 

-k  -k 

Then  we  have  x_  = A.  x,  = A,  x.  , - ™ < k < <»  which  yields 

1 K1  ’ *2 


(b)  The  condition  c follows  as  in  (a).  The  converse  can  be 
proved  by  constructing  y^ . Multiplying  both  sides  of  equation  (2. 4)  by 
C,  one  obtains 

Cxj.  = CA^kx  r C (fj)  since  A^x  r f1  (2.8) 
but  comparing  equations  (2.t)  and  (2.8),  one  obtains  that 
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1 

f 


! I 

I T 


5 


7 


I 


yk  = CA^y  £ C(f1 ) . i 

Uniqueness  follows  from  the  uniqueness  of  x. 

The  following  examples  illustrate  Lemma  2: 

i 

Example  2.1 : 

Consider  homogeneous  LSC  defined  over  G^(2)  characterized  by  the 
following  set  of  equations. 


— —1 
X1  ( k+1 ) 

‘ 1 

1 

— - 

xlk 

X2( k+1 ) 

L J 

0 

1 

x2k 

L. 

>'k 

■ i 

ni 

’] 

1 

'*u1 

(2.9) 


(2.10) 


with  *o  ■*•[*]«  fl  • 

■[:]■ 


Clearly  xk  = 


x„. 

r k 


f ] Vk  and  yk  = 1 V*  i.e.,  yk  c C ( f n ) 


Now 


consider  xq  = x - ^ j . f(,  which  yields. 


(xk)  = { 


0 10  10  1 
’ 1 • 1 ’ 1 • 1 ’ 1 ’ 1 


1.  i.e.. 


xk  c f ^ for  all  k and  the  output  sequence,  lyk),  with  y - 1 , is 

f yk i * i , l,  o,  i,  o,  l,  o, } 

i.e.,  yk  c C ( f 1 ) for  all  k. 

in  the  above  example,  note  that  the  uniqueness  of  {yk">  is  solely 
determined  by  the  "central  value"  of  the  state,  xq  = x_. 

In  the  following  example,  two  distinct  central  values,  xo  = x, 


A 
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one  in  f^  and  the  other  not  in  fj  are  considered.  It  is  then  shown  that 
the  solution  exists  for  x c f^  but  does  not  exist  for  x not  in  f ^ . 
Example  2.2: 

Consider  a homogeneous  LSC  defined  over  GF(?)  characterized  by 


X1  (k+1 ) 

~1 

1 ' 

xn 

*2(k+l)  _ 

0 

0 

- 

x2k 

- 

Clearly, 


Consider 


R(A)  = R(A2)  * R(A3)  = = fv 


Then  A1  : f1  f ^ where  A^ 


With  central  value 


*■[!]• 


* t i] 


one  obtains, 


for  al 1 k 


Now  let  x = |j  j notjn  f.j . Then  one  obtains 


■[*] 


not  in  f ^ . 


Lemma  3:  (Non  Homooeneou’  Case) 

Consider  the  t-quations,  xk  + ^ = Ax^  + Bu^,  - 00  < k < «*•  , defined 
over  a finite  field  and  let  u^  be  a periodic  sequence.  Then  there 
exists  x such  that  the  equation  has  a periodic  solution,  (x^i, 

x " x . 

0 


Proof : 


It  is  well  known  that  when  a finite  state  machine  is  driven  by  a 
one-sided  ultimately  periodic  input  sequence,  the  output  sequence  (and 


therefore  state  sequence)  produced  by  the  machine  exists  for  any  ini- 
tial condition  and  is  ultimately  periodic.^ 

Let  { uk ) , k = 0,  1 , 2 , . , . . be  the  "one-sided"  ultimately 
periodic  input  sequence  constructed  from  a two  sided  periodic  sequence 
{uR}  via  uR  - uR , k 0.  Let  (xR)  ~ fx^,  * ^2*  ♦ • • • * x^-j-t 
. . . . , xMT_1 , . . . (where  xq  is  any  initial  state)  be  the 
resultant  "one-sided"  ultimately  periodic  state  sequence  produced  by 
«u.  ■ . Without  loss  of  generality  one  may  assume  that  the  values  x^j, 


'NT+1 


iMT.r  M 'N,  constitute  one  period  of  the  periodic 


portion  of  { x . where  T is  the  period  of  'uR'».  Next  define  a periodic 
sequence  (y.  * which  coincides  with  ix  ■ for  NT  * k < MT-1  i.e.,  xR  = 
where  r is  the  unique  integer  r * k ♦ t;(M-N)T  and  NT  <_  r *_MT-1,  q is 
an  integer  (positive,  zero  or  negative''.  !t  may  be  shown  that  the 
sequence  •.  xR  • satisfies  t nonhoio  igneous  equation. 

Case  1 : (NT  < r V'-?'. 

Substituting  r - k * in  the  nutation  * ^ = Ax^.  + Bu^  yields 

Xk  + l+q(MT-NT)  = AxJ *q(V’-\T)  + Puk4q(MT-fjT) 


. o'"’  -*,T  1 + 


Which  is 


Xk  + 1 A\  + S 


Case  2-  (r  - MT-1 ) 

Substituting  r - MT-1  in  the  equation  x .j  = Ax^  + Bur  yields. 


XMT  = AxMT-1  + BuMT-1  ’ S1nrC  uMy-l  = UNT-1 


one  obtains , 
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or  x^  3 Ax^  + Bu^  for  k = NT-1. 

Finally  note  the  uniqueness  of  a periodic  portion  of  {x^l  i.e.t 

XNT*  XNT+1  * ’VlT-l  where  T is  the  Period  of  tufc},  implies 

the  uniqueness  of  {x^}  satisfying  x^+^  = Ax^  + Bu^,  - « < k < «. 

The  following  example  illustrates  Lenina  3: 

Example  2.3: 

2 

Consider  an  ISC  defined  over  GF(2  ) and  characterized  by  the  set 
of  equations, 


(2.1?) 


- 

xl(k+l) 

1 

1 

X 1 ( k ) 

— ■» 
0 

= 

X2(k+1 ) 

0 

0 

_x?(k)  ^ 

1 

Let  {ukl  = {1,  a,  l+a.  a,  1 +a , 


and , 

-fil 


which  give 


'V  ‘ |0 


Then 


X1 

x2  *3 

x4 

x5 

xe 

0 

1 l+a 

0 

1 

1 + C1 

1 

a 1 +a 

1 

a 

1 

periodic  portion 

of 

{xk 

> such 

tnat 

x3  = XNT 

1 

1 

l+a 

0 

1 

l+a 

t 

* 

i 

> 

* 

* ♦ 

o 

a 

l + i 

1 

n 

1 +a 

where  x 


• ' [::] 


is  the  solution  to  equation  (2.1?). 


2.2  Spectral  Theory: 

For  an  operator  H on  a vector  space,  one  says  that  a scalar,  >,  is 
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said  to  be  an  eigenvalue  if  there  exists  a non-zero  vector  l)  such  that 
MU  = M (2.13) 

and  one  says  that  U is  the  corresponding  eigenvector.  Although  this 
definition  is  traditionally  associated  with  real  or  complex  matrices 
on  a finite-dimensional  space,  the  defining  equality  for  eigenvectors 
and  eigenvalues  also  holds  for  arbitrary  operator  on  an  abstract 
vector  space.  In  particular,  one  can  use  the  above  equality  to  define 
the  eigenvalues  and  eigenvectors  of  a single  input  - single  output 
LSC  by  viewing  it  as  a linear  operator  on  the  infinite  dimensional 
vector  space  of  "two-sided"  sequences  with  values  in  GF(p).  In  this 
the  defining  equality  becomes. 


Vi  = Axk  4 Buk 


;uk  r Cxk  4 Juk 


(2.14) 


Here  {u^'-  is  an  etgenseauence  takinn  its  values  in  a finite  field 
and  the  eigen  values  in  general  , takes  their  value  in  the  algebraic 
closure  of  GF(o).  This  is  easily  verified  by  observing  that  the  eigen 
values  for  a matrix,  M,  over  GF(p)  are  the  zeros  of  che  polynomial 
det(-I-M).  Althoug'  the  coefficients  of  this  polynomial  are  GF(p), 
these  zeros  of  the  polynmna  may  lie  in  its  algebraic  closure.  The 
algebraic  closure  of  a finite  field  plays  essentially  the  same  role  for 
an  LSC  defined  over  that  field  as  complex  numbers  do  for  an  analog 
circuit  defined  over  the  real  field.  The  minor  problems  arising  due 
to  the  use  of  the  algebraic  closure  are  discussed  in  the  next  section. 

For  th«  periodic  input  = iek;,  k = 0,  si,  ±2,  ....  (2.15) 

where  e is  an  element  of  the  algebraic  closure  nf  the  finite  field 
over  which  an  ' Sf  is  defined.  Cefine  a "transfer  function"  in  the 
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usual  manner, 

S(e)  - .1  + C (el-A)'1  B (2.16) 

Here,  S(e)  is  a rational  function  in  e with  coefficients  in  a 
finite  field,  whose  poles  and  zeros  are  well  defined  elements  of  the 
algebraic  closure  of  the  finite  field.  For  each  e (except  for  the 
poles  of  S(e)),  the  sequence  (u^l  * (eS.  - - < k < *>  is  an  eigen 
sequence  of  the  LSC  with  eigenvalue  S(e).  The  eigenvalues  of  the  LSC 
are  elements  of  the  algebraic  closure  and  the  elgensequences  take  their 
values  in  the  algebraic  closure  just  as  eigenvalues  and  eigenfunctions, 
for  a system  defined  over  the  real  field  are  complex.  The  fundamentally 
important  fact  here  is  that  the  eigensequence  is  independent  of  the 
LSC  under  study  and  hence  is  assured  to  be  an  eigensequence  for  a 
faulty  circuit  even  if  the  fault  is  unknown. 

Theorem  1:  Let  an  LSC  be  characterized  by 

■Vi  • Axk  * Buk  (2-1A) 

- OD  < k < OO 

yk  = Cxk  + Juk  (2- IB) 

over  a finite  field.  Then  for  each  e in  the  algebraic  closure  of  the 
finite  field  for  which  S(e)  = J + C (le-A)"1  B is  defined,  there 

exists  an  initial  condition  x for  (2.1)  such  that  the  sequence 

c 

k 

(u^l  = ie  j(  k = 0,  ±1,  + 2 , is  an  eigensequence  f°r  LSC 

with  this  initial  condition  and  eigenvalue  is  S(e). 

Proof: 

4 lr 

Since  { u^}  = {e  },  k = 0,  -1,  *2 is  a periodic  sequence, 

the  existence  of  a solution  to  (2.1)  is  guaranteed  by  Lemma  3.  The 

k 

existence  of  an  initial  condition  such  that  -f 1 = {e  } , k = 0,  t.l , *2,  . 
. . is  an  eigensequence  with  the  eigenvalue  3(e)  can  be  proven  by 


substitution  as  follows: 
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The  claim  is  that  - (Ie-A)'^  Be*'  (2.17) 

is  a unique  solution  to  xk+^  = Ax^  Bu^,  - « < k < « (2.18) 

k 

with  an  input  sequence  i uk } = {e  },  k = 0,  ±1 , i2,  .....  such  that 


*0  = ( Ie-A) B.  To  show  this,  substitute  (2.17)  in  (2.18) 
AxR  + Buk  = A[(Ie-A)_1  Bek]  + Bek 
- ( Ie-A)-1  Bek  ♦ 1 
= xk+1  ' 


Uniqueness  of  xk  can  be  shown  be  contradiction,  let  there  be  two  solu- 
tions x. ' and  x."  to  equation  (2.18)  such  that  x ' - x " = (Ie-A)~^  B. 
k k oo 

Then  one  obtains 


vr  r av  - Buk 

(2.19) 

Vl"  ‘ AXk"  * BUjc 

(2.20) 

Subtracting  (2.20)  from  (2.19) 

<vr  - = A<v  - v> 

(2.21) 

Let  x(  ‘ - xk"  = xk  then  (2.21)  becomes. 


Then  from  l emu  one  nets 


which  implies. 


- 1 k 

This  shows  xk  = (Ie-A)  Be  is  the  unique  solution  to  equation  ( 2 . 1 A ) 
such  that  xQ  = x = (Ie-A)-1  B. 

Substituting  (2.17)  into  (2.16)  yields  • 


»k = c\  * Juk 
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= [C(Ie-A)'1  B + J]  ek 

= S(e)  ek  (2.22) 

Uniqueness  of  y^  follows  from  the  uniqueness  of  x^  \tk. 

The  essence  of  Theorem  1 is  that  it  allows  one  to  interpret  the 
"transfer  function"  of  an  LSC  as  a function  such  as  is  done  for  contin- 
uous time  systems  over  the  real  field  rather  than  as  an  abstract  oper- 
ator as  is  usually  done  for  LSCs.  Thus,  it  is  a function  which  iden- 
tifies an  eigenvalue  with  its  eigenvector.  In  the  case  of  multiple 
input  -multiple  output,  the  above  arguments  go  *l.~ough  with  the  matrix 
S(e)  = J + C(Ie-A)'1  B interpreted  as  a matrix  of  eigenvalues. 

Example  2.4: 


f i r\n  1 i 


innnf-cinnlfi  ftiitnut  I inOAr  ^AniiPfit  T A 1 fir- 

» ll^w  V • I • J « WW1  w “ . . • w — 1 • 


cuit  shown  in  Figure  2.2  and  described  by  the  equations 
xk+l  = xk  + uk 


(2.23) 


defined  over  GF(2). 


Substituting  the  values  of  A,  B,  C,  J,  into  (2.16)  yields 

Me)  = y-h- 

Consider  the  extension  field  GF (2 } whose  elements  are  (0,  1,  a. 


(2-24) 


23 


S 


4* 


i 


l 


1 


! 


r 


» 


a + 1 = a2)  and  let  e = a,  then  k.  { u k } - {ek}  = {ak},  i.e...  uQ  = 1 

and  { Uk)  = {.  . . . 1 , ft,  1 + a,  1 , a,  1 + o,  1 t a,  1 + c,  . . . .} 

k k 

which  produces  an  output  sequence,  = S(e)e  - S(a)a  „ i.e.,  ly^l  s 

{.  . . . a,  a2,  1 , a,  a2,  1 , a,  a2,  1 ) where  yQ  = a . Now 

L 

letting  e = a + 1,  {uk>  = {a  + 1}  produces  an  output  sequence  (y^i  = 

{.  . . . a2,  a,  1,  Cl2,  a,  1,  a2,  a,  1,  . . . .)  where  yQ  = a2. 

k 

In  the  above  example,  the  eigensequence  {u^}  - {e  1 is  sso- 
ciated  with  an  eigenvalue  S ( e ) = p — . This  example  sho'-'S  that  by 
choosing  different  elements  from  the  extension  field  of  6F(2),  differ- 
ent output  sequences  are  generated.  Also  note  that  in  this  case  the 
state  and  the  output  sequences  will  be  the  same. 

2.3  Computing  S(e)  from  the  Impulse  Response: 

Using  input  signals  taking  their  values  in  the  algebraic  closure 
of  a finite  field  presents  two  problems.  First,  one  must  be  able  to 
do  computation  in  the  extension  field  and,  second,  one  must  actually 
test  the  system  with  signals  that  take  their  values  in  the  original 
field  since  t.he  physical  system  is  not  capable  of  accepting  inputs  from 
the  extension  field.  The  first  problem  can  be  solved  without  much 
difficulty  since  GF(pn)  car,  be  represented  as  a field  of  rntkl  degree 
polynomials  with  coefficients  in  GF(p).  Here,  addition  is  the  usual 
polynomial  addition  and  multiplication  is  the  usual  polynomial  multi- 
plication modulo  of  the  (m  + l)s^  order  irreducible  polynomial. 

The  second  problem  can  be  solved  by  simply  measuring  the  impulse 
response  of  the  system  and  then  computing  the  "exponential"  response. 

Since  the  goal  here  is  to  be  able  to  find  the  relationship  be- 
tween the  zero  state  impulse  response;  { h k 1 , k >_  0,  and  the  matrix 
transfer  function  S(e),  one  can  assume  that  the  input-output  relation- 
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J 


I 

t 

I 


1 1 


2‘ 

ship  of  the  LSC  is  characterized  by  a minimal  (A,  B,  C „ J)  realization, 
which  characterizes  the  actual  physic?!  system  under  study.  Since  such 
a realization  is  controllable  it  can  always  be  driven  from  the  zero 
state  at  time  k = kQ  to  any  desired  state  at  time  k = 0 by  a sequence 
of  inputs  b_j,  kQ  j < 0 provided  kQ  <_  - 6 where  6 is  the  degree  of 
the  LSC. 

The  relationship  between  the  "exponential"  response  of  a single 
input  - single  output  LSC  and  its  zero  state  impulse  response  is 
given  by  the  following  theorem.  The  results  of  this  theorem  can  easily 
be  generalized  for  the  case  of  multiple  input  - multiple  output  LSC's. 
Theorem  2:  Let  k ^ 0 be  the  zero  state  impulse  response  of  a 

single  input  - single  output  LSC  and  S(e)  = J + C(le-A)"^  B be  its 

transfer  function.  The  "exponential  response"  $(e)e  to  an  input 

, . k.  . 

sequence  i > = ie  i is  yiven  by, 

k k°  k'1  k i-i 

S(e)ek  = l h,  , b . + l h.  , ek  1 , k > 0 (2.25) 

j=l  J k i=-l  11 


The  b .'s  are  (unknown)  inputs  that  drive  the  system  from  the 
~ J 

zero  state  at  time  k = kQ  to  xQ  - x (possibly  non-zero)  at  k = 0. 
Proof: 


(e 


k < -1 


Define  { ek } - l 


(2.26) 


(e  k)  $ 


k > -1 
k < -1 


(2.27) 


k > -1 
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Then  { ek>  * ieKl  + {e  k} , k = 0,  ±1 , tZ Let  yfc  be  the 

response  of  the  LSC  to  {u^}  = {e  1.  Let  ^ ^ be  the  response  of  the  LSC 

|r  _ 

due  to  {u^  = {e^  }.  Let  be  the  response  of  the  LSC  due  to  {u^}  = 
k, 

ie  } . 

It  follows  firectly  from  the  defining  equations  for  the  minimal 
realization  of  the  LSC  that 


, k-1  . i.  - l 

x.  = Akx  + ][  A1  B e 1 

K 1=0 


(2.29) 


and 


k-1 


y = CA*x  + £ CA  B e 

* i =0 


i d „k-i-l  + j ek 


(2.30) 


0 be  a sequence  of  inputs  that  urives  the 


Mow  let  b_j , kQ  •-  -j 
minimal  realization  from  the  zero  state  at  time  k = to  x at  time 
k = 0.  Such  a sequence  exists  for  a minimal  realization  if  kQ  < -6, 
where  5 is  the  dimension  of  the  minimal  realization. 

It  follows  that 


'-k 


= 0 


x-(k  -1)  = Bb-n 


x-(ko-2)  ” ABb-n  * Bb-(n-l ) 

Mk  -3)  ‘ A?Bb-„  * + Bb-(n-2) 


I-  *o‘  x-(k  -k  ) 

' O 0 


= [B  AB  . . 


„ k -?n  .k  - 1 -j 

A o B A o BJ 


-2 


J-(k  ) 
' o' 


V, 
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I! 


which  gives 

CAkx  = [CAkB  CAk+1B 


CAk+ko“2B  CAk+ko’]B] 


’-1 


b-(k0-D 

b-k 


(2.31) 


The  values  of  the  Impulse  response  {h^}  are  given  In  terms  of  the  ISC 
minimal  realization  description  by. 


h0  = 0 
h]  = C8 
hz  = CAB 


CAk_1B 


k > 0 


Substituting  (2.32)  In  (2.31)  yields, 
CAkx  - [hk+1  hR+2  ... 


hk+k  -2  hk+k  -1^ 
O 0 


J-1 


-2 


(ko-D 


o 


(2.32) 


‘ il  hj*k  b'-« 

Substituting  (2.33)  into  (2.30)  yields, 
k k] 

yk  - 1°  hj.k  t>_i  + I CA1  B e 

K j-i  3 11  J i=l 


k-i-1  . , k 
+ j e 


(2.33) 


Utilizing  (2.32)  yields. 
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ko 

k-1 

I,  Vk  b-j + 

l 

i=0 

k<, 

k-1 

I,  Vk  h-J  * 

y 

i=-1 

From  (?  2)  the  "exponential  response" 

yk  - S(e)ek 

Thus  from  (2.34).  it  follows  that, 


k-i-1  . k 

e + hQe 


h1+1  e 


k-i-1 


is  given  by 


(2.34) 


S(e)ek 


k-1 


l 

i=-l 


h . , , e 
i+l 


k-i-1 


In  ; s where  the  ISC  is  not  in  the  zero  state,  the  zero  state 

impulse  rev  se  can  be  determined  as  follows:  Observe  the  response 
of  the  LSC  due  to  an  input  sequence  lu^1  = {01.  When  this  response 
becomes  periodi  , apply  an  impulse  sequence  to  the  LSC.  The  zero 
state  impulse  «pon$e  is  then  given  by  the  dii  Terence  of  two  responses 
i.e.,  responses  obtained  after  and  before  the  impulse  sequence  is 
applied. 

To  apply  Theorem  2,  one  observes  that  Equation  (2.25)  is  linear 

in  all  the  unknowns  (S(e)  and  the  b_j,  ko  -j  <_  0)  hence  by  writing 

the  equation  for  1=1,2,....,  ky+l , one  may  set  up  a kQ+l  by 

kQ+l  matrix  equation  which  may  be  solved  for  S(e)  in  terms  of  the 

measured  values  of  h.  , 0 < k < 2k  +1.  In  the  case  of  multiple  input  - 

k o 

multiple  output  LSC's,  S(e)  is  a matrix  of  eigenvalues  and  the  impulse 
response  is  also  a matrix.  The  procedure  for  determining  S(e)  in  the 
single  input  - single  output  case  can  easily  be  generalized  for  deter- 
mining the  impulse  response  matrix  and  (in  turn)  the  matrix  S(e). 

The  following  example  illustrates  Theorem  1. 
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Example  2,5: 

Consider  the  ISC  defined  In  Example  2.3.  The  zero  state  Impulse 
response  Is  given  by 

h0  * 0,  h1  = h?  = h3  =• - 1, 

from  Theorem  1 

k 0 k- 1 k i 1 

S(e)e  - ^ hJt,  b_j  + J ^ h,+,  e ' ' 


Let  e = a e GF(22),  k = 1. 

For  k = 1,  S(a)a  = b_1  + 1 (2.35) 

and  for  k * 2,  S(a)a2  = b_^  + a + 1 (2.36) 

Solve  (2.35)  and  (2.36)  simultaneously,  one  obtains, 

S(«)  = a - y-]—  (2.3?) 

o ? 

Now  let  e = o'  - 1 + a e w(2  )* 

for  k * 1,  Sta2)^2)1  = + 1 (2.38) 

for  k * 2,  S(a2)(a2)2  « b}  + a2  + 1 (2.39) 

Solving  (2.38)  and  (2.39)  simultaneously,  one  obtains, 

S(a2)  = 5-  (2.40) 

1 + a 


Equations,  (2.37)  and  (2.40),  could  be  verified  by  substituting  the 
corresponding  value  of  e in  equation  (2.24). 


CHAPTER  III 


FAULT  ANALYSIS  - LINEAR  SEQUENTIAL  CIRCUITS 
3.1  Component  Connection  Model: 

A Linear  Sequential  Circuit  Is  usually  characterized  by  ar.  Input- 
output  state  model  for  the  purpose  of  fault  analysis;  wherein  one 
attempts  to  collect  Input-output  data  to  determine  faulty  components 
within  the  LSC.  The  fault  analysis  algorithm  developed  in  this  chapter 
use*-  the  component  connection  which  relates  input-output  behavior 
directly  to  component  parameters  rather  than  the  state. 

The  component  connection  model  was  first  intuitively  used  by 
1 8 

Prasad  and  Traboth  and  has  been  used  by  several  other  investigators 
in  the  area  of  fault  analysis  in  analog  circuits.*’5,14 

The  primary  reason  in  choosing  the  component  connection  model  for 
fault  analysis  is  that  it  is  so  heavily  algebraic  that  it  unities 
various  graphical  and  diagramatic  connection  theories  and  at  the  same 
time  smooths  the  transition  from  mathematical  model  to  computer  algo- 
rithm. 

In  the  sequel  the  digital  version  of  the  component  connection 
model  is  developed  which  is  essentially  the  same  as  in  the  analog  case 
except  for  minor  changes.  The  development  is  repeated  here  for  the 
sake  of  completeness  and  for  demonstrating  the  interpretation  of  terms 
used  in  the  model  for  the  digital  case. 

A mathematical  interpretation  of  a system  is:  a mapping  from  a 

set  of  inputs  to  a set  of  outputs  i.e.,  an  input-output  relation. 
Letting  u and  y represent  system  input  and  output  sequences,  respec- 
tively, with  values  in  a finite  field  or  its  algebraic  closure,  we  may 


29 


abstractly  denote  a system  by 
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y = Su 


(3.1) 


where  S is  the  mapping  (operator)  representing  the  system.  The  system. 
S,  is  primarily  determined  by  two  factors:  the  component  types,  and 

the  ways  the  components  are  interconnected.  In  the  case  of  an  ISC,  a 
system  consists  of  a fixed  interconnection  of  several  Linear  Sequential 
Circuits,  each  in  turn  which  is  a fixed  interconnection  of  components 
or  devices.  Components  can  be  scalars,  adders,  delayors  or  linear 
gates.  The  form  of  the  algebraic  connection  model,  known  as  the  Com- 
ponent Connection  Model,  can  be  conceptually  deduced  f om  Figure  3.1. 


Figure  3.1:  System  Representation 

In  Figure  3.1,  u and  y represent  overall  system  inputs  and  outputs 
(finite  field  valued  sequences)  while  a and  b are,  respectively,  com- 
ponent input  and  output  sequences  over  the  same  field. 

By  a "Component  Connection  Model"  is  meant  a system  model  where- 
in the  components  and  connections  are  characterized  by  separate  equa- 
tions. In  particular,  one  assumes  that  each  component  in  a system  is 


characterized  by  the  D-transform  equation. 


b-(D)  = [ (g. )f i (D)]  a • (D) 


i = 1,  2,  . . . n (3.2) 


■i 
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where  a^D)  and  b,(D)  are  respectively  the  D-transtforms  of  the  Input 
and  output  sequences  of  the  I**1  component.  The  dynamics  of  the  1 
component  are  represented  by  f ^ (D ) and  the  scalar  constant,  9j» 
gain  of  the  I**1  component.  For  the  purpose  of  fault  analysis.  It  Is 
assumed  that  the  component  dynamics,  f^{D)v  do  not  change  even  after  a 
fault  has  occurred  and  all  failures  manifest  themselves  as  changes  in 
with  memory  elements  and  connections  good.  This  guarantees  that  a 
linear  system  falls  linearly.  Such  a fault  model  includes  the  usual 
"open"  and  "short  circuit"  faults  and  "stuck-on-zero"  faults.  "Stuck- 
on-one"  faults  can  not  be  Included  in  a strictly  linear  theory  since  a 
device  which  is  "stuck-on-one"  is  nonlinear.  "Stuck-on-one"  faults  are, 
however,  included  in  the  generalization  of  the  theory  to  Affine  Sequen- 
tial Circuits  described  in  Chapter  V. 

In  actual  practice  one  normally  works  with  the  n separate  component 
equations  like  (3.2).  hotationally  this  may  be  combined  into  the 
single  matrix  equation, 

b(D)  = [(G)F(D)]  a(D)  (3.3) 

where  b(D)  = Col.  ( ( D ) ) 

a(D)  = Col-  (a^D)) 

G = diag.  (G.) 

F(D)  = diag.  (^(D)) 

To  obtain  a mathematical  model  for  the  connections,  redraw  Fig. 

3.1  as  in  Fig.  3.2  where  the  components  and  connections  are  shown 
separately.  The  connections  may  be  viewed  as  a separate  multiple 
input  - multiple  output  component. 


32 


Components 


rn 


Connections 


Figure  3.2:  Component  Connection  Model  Representation 
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Since  connection  elements  are  linear  and  algebraic,  and  connections 
can  be  characterized  by  linear  algebraic  constraints  (Adders,  scalers, 
etc.)  the  connection  model  car*  be  represented  by  the  matrix  equation. 


(3.4) 


Entries  in  the  l matrices  take  their  values  in  the  same  finite  field 
over  which  the  ISC  under  study  is  defined.  For  digital  circuits  the 
t..  matrices  are  usually  permutation  matrices  describing  how  the  out- 
puts  of  one  component  are  connected  to  the  inputs  of  another.  In 
Equation  (3.4),  u(D)  and  y(D)  ere  respectively  the  D-transforms  of  the 
externally  accessible  inputs  and  outputs  of  the  LSC.  The  digital  ver- 
sion of  the  component  connection  model  is  sufficiently  general  to  in- 
clude most  ISC's,  although  it  is  not  universal-  Equation  (3.4)  has  the 
symbolic  representation  shown  in  Figure  3.3. 


Components 


Figure  3.3:  Symbolic  Representation  of  Matrices 

In  the  following  example,  L.^,  G and  F(D)  matrices  are  determined 
for  the  LSC  shown  in  Fig.  3.4. 

Example  3.1 : 

Consider  the  system  defined  over  GF (2)  as  shown  in  Figure  3.4. 


Figure  3,4:  Figure  for  Example  3.1 


The  L matrices  for  the  system  described  in  Figure  3.5  are  given  by 


T-»*> 


3' 


Simultaneous  solution  of  equations  (3.1),  (3.3)  and  (3.4)  yields  a 
complete  description  of  the  LSC  input-output  relationship  of  the  form, 
y(D)  = [5(D)]  u(D)  (3.5) 

= U L 22  + L21(l-(G)F(°)Li1)-1(G)F(D)L12]u(D)  (3.6) 

That  is,  the  overall  system  operator  is  given  in  terms  of  the  com- 
ponent dynamics,  F(D),  component  parameters,  G,  and  connections,  l-., 

* J 

by  the  equality, 

5(D)  L22  + L?1(1-(G;F(0)L11)'1(G)F(D)L12  (3.7) 

S(D)  is  a matrix  of  rational  functions  in  D with  coefficients  in 
a finite  field,  whose  poles  and  zeros  are  well  defined  elements  of  the 
algebraic  closure  of  the  finite  field. 

In  view  of  the  hypothesis  that  all  t SC  failures  manifest  them'. elves 
as  changes  in  G,  with  F (D ) and  L^.  remaining  constant,  it  is  natural  to 
view  Equation  (3.7)  as, 


s 


t. 


*<* 


c 


l- 


? 


i 

I Jt 


I 
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S(D)  = fD(G)  (3.8) 

where  fp  is  a nonlinear  function  even  though  S Is  a linear  operator. 

The  function,  fp,  Is  entirely  determined  by  the  component  dynamics  and 
the  connections  relating  component  parameter  values  to  the  system's 

5 

input  and  output.  The  function,  fp,  is  called  the  connection  function. 
Equation  (3.8)  is  in  just  the  right  form  for  the  study  of  fault  analysis 
because  faults  are  assumed  to  manifest  themselves  only  as  changes  in  G. 
3.2  Faul t Analysis: 

The  following  theorem  gives  the  relationship  between  S(e)  which 
may  be  computed  from  the  zero  state  impulse  response  of  the  LSC  and 
$(D). 

Theorem  3: 

Given  an  LSC.  described  by  y(D)  = S(D)u(D).  Let  { u^,}  = {e+k},  k = 0, 
±i,±2,  . . . . , where  e is  an  element  in  the  algebraic  closure  of  a 
finite  field  over  which  the  LSC  is  defined.  Then  S(e)  j D = e (3.9) 
Proof: 

Consider,  first,  the  single  input  - single  output  case.  In  this 

case  S(D)  is  a rational  function  in  D with  well  defined  poles  and  zeros 

that  can  be  expressed  as 

a + a.D  + a0D"  + — + — + a D 
ol2  n 

S(D)  (3.10) 

1 + b.D  + bo0^  + — + — + b Dm 
1 2 m 

which  implies 

a + a.D  + a„D2  + + a Dn 

ol2  n 

y(D)  = u(D) 

t + b.O  + b,D2  + + b/ 

1 2 m 


Equivalently, 
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y(D)[l  + b^D  + b202  + — + bmOm]  = [&Q  + + a2D2  + — + a^^uCD). 

(3.10A) 

Since  Dmy(D)  - {yk+Tn}  and  Dnu(D)  = {u^},  (3.10A)  becomes, 

{V  + VW  + — + VW  = ao{uk}  + Wl>  + ~ + VW 

(3.11) 

Let  uk  - {ekl  , k = 0,  ±1,  ±2,  .... 

Then  yk  = {ye+k}  (3.12) 

Substituting  (3.12)  into  (3.11)  yield, 

{ ye k } [ 1 + b^e  + b2e2  + — + b[r)em]  = {ek}[aQ  + + — + anen] 

or  a0  + a^e  + a2e2  + — + a/ 

{UjJ 

1 + b,e  + b^e2  + — + b em 
1 2 m 

= S(e){uk>,  where  S(e)  = S(D)  | D = e (3.13) 

The  multiple  input  - multiple  output  result  follows  easily  by  applying 
the  above  analysis  over  all  possible  input-output  pairs  with  all  inputs 
zero  except  the  one  corresponding  to  the  pair  in  question. 

Given  the  Theorem  3,  Equation  (3.6)  may  now  be  written  as 
S(e)  - ffi(G)  = [l2?  + L21d-(G)F(e)L11)'1(G)F(e)L12].  (3.14) 

If  (3.14)  can  be  solved  for  G,  the  fault  analysis  of  the  Linear 
Sequential  Circuit  is  complete.  Unfortunately,  however,  most  LSC's 
have  more  components  than  input-output  equation,  in  which  case  (3.14) 
represents  a set  of  equations  that  has  more  unknowns  than  the  number 
of  equations.  However,  this  difficulty  can  be  overcome  by  exploiting 
the  LSC  dynamics.  Indeed,  this  exploitation  of  the  circuit  dynamics  is 
the  key  to  the  spectral  theoretic  approach  because  it  permits  the 


number  of  externally  accessible  LSC  terminals  required  for  fault  analy- 
sis of  an  LSC  to  be  reduced  from  that  required  for  a combinatorial  cir- 
cuit of  similar  complexity  rather  than  increased  as  Is  the  case  with 
traditional  fault  analysis  techniques.  To  see  this,  note  that  Equa- 
tion (3.14)  is  valid  for  any  e in  the  algebraic  closure  of  a finite 
field  over  which  the  LSC  is  defined  and  that  <he  resulting  equations 

are  dependent  on  the  choice  of  e.  As  such,  more  equations  can  be  creat- 

ed without  changing  the  number  of  unknowns  by  writing  the  set  of  simul- 
taneous equations  for  different  e's  in  the  algebraic  closure  as, 

S(e1 ) = f (G)  = [L22  ^ L21(l-(G)F(e1)L11)'1(G)F(e1)L12] 

S(e2)  = fe  (G)  = [L22  + L21(l-(G)F(e2)L11)"1(G)F(e2)L12] 

I 2 (3.15) 

s(ek)  = f (G)  - [L22  + L21(l-(G)F(ek)L11)"1(G)F(ek)L12]. 

In  the  case  of  a combinatorial  circuit  f is  independent  of  e. 

ei  1 

hence  the  additional  equations  are  not  independent  and  do  not  simplify 
the  fault  analysis  process.  In  the  case  of  a Linear  Sequential  Circuit, 
Equations  (3.15)  will  be  independent  and  solvable  for  G even  though  a 
single  Equation  (3.14)  is  not  solvable. 

With  the  above  theory,  one  can  bormulate  the  fault  analysis  algo- 
rithm. The  procedure  consists  of  the  following  steps. 

(i)  Measure  the  zero  state  impulse  response  of  an  LSC  under  study. 

(ii)  Compute  S(e^)  from  (2.25)  for  various  elements  e-| , e2,  e^,..., 

in  the  extension  of  the  finite  field  over  which  the  LSC  under  study  is 
defined.  The  number  of  e's  that  one  should  choose  depends  on  how  many 
component  parameters  are  to  be  solved. 


38 


* 

(ili)  From  $(ej),  S(e2) S(efc)  calculated  In  Step  (11)  and  from 

the  known  matrices  L.,,  write  the  simultaneous  equations  (3, IS). 

■ J 

(1v)  Solve  the  equations  (3.15)  for  6, 

In  the  next  section,  the  algorithm  for  solving  (3.15)  Is  formulated. 

3.3  Equation  Solving: 

A procedure  is  now  given  for  solving  the  set  of  simultaneous 

equations  (3.15).  This  procedure  is  based  on  a "term  expansion"  algo- 

14 

rithm  developed  by  Ransom  for  analog  fault  analysis, 
let  T(D)  - {(G)F(D))*\  Equation  (3.7)  becomes, 

S(D)  = L22  + L21(T(D)-L1,)"1L12,  i = 1 ,2, . . . ,k  (3.16) 

Next,  perform  the  "term  expansion"  of  the  "inverse"  in  (3.16). 

(T(D)-L^)  = — adj.  (T(D)-L^ ) where  t denotes  transposition,  A = 

(T(D)-L^),  and  adj.(T(D)-L.^)  denotes  the  matrix  whose  i-jj^  element 

is  the  ijjth  cofactor  of  (T(D)-L^). 

Taking  che  "vec"  operation  on  both  sides  of  (3.16)  yields 

vec(S)  = vec(l2?)  + (L^*  ® <-21)  v~  GaPa  (3.17) 

3 3 

19 

where  © denotes  the  Kronecker  matrix  product. 

In  arriving  at  (3.17)  one  needs  to  use  the  formula  vec(ABC)  = 

(CT  © A)vec(B),  the  equalities  a = Vo  and  vec  adj. (T-L.. )*  = Go 

3d  M da 

where  V is  a row  vector,  G is  a matrix  and 


Reproduced  from  3( 

bett  available  copy.  wdP 

p ■ ■ i i — 

4 

Here  n is  the  dimension  of  T(D}  and  T is  a vector  whose  elements 

are  all  product  combinations  of  the  diagonal  elements  of  T taken  i at  a 

time.  That  is,  for  an  LSC  with  three  g's,  = [1,  , c^,  g^,  g^, 

g,g0,  q0g0,  g,q0g0].  Here  the  elements  of  V.  and  G are  constants  in 

an  extension  of  a finite  field.  An  alcoritbm  for  determining  V,  and  G 

a d 

is  given  in  reference  and  will  not  be  repeated. 

One  can  write,  V r>  = Vp  + r 

d a 

Gc=G,jt).  (3.18) 

where  p,  V,  G are  o.  , V and  G , respectively,  with  their  first  element 

d S.  cl 

(column)  deleted.  The  first  t lenient  (column)  of  V,  and  G are  respec- 

a n 

lively,  r and  . 

Substituting  (3.1^)  into  (3.17)  yields. 


■:(lr*  © l91)g  - vec(.T-:.„;v;  £.  = -'l.,.*  ® u,)>. 


'ec (T-1^2 ) r 


Writing  Equation  (3.13)  i y '-'ariour,  vduc^  of  e- , one  can  obtain 
the  desir'd  form  of  set  if  r i:.:u ; Laneui.  '•>  equations, 

«•■>,)  i = ‘(f,) 

i • 

*(e,)  . - Me..) 


C ( e k ) >•  r '(<%.' 

wr.ich  can  be  written  mere  compactly  with  the  obvious  notation?!  defin- 
itions as. 


In  general,  may  have  linearly  dependent  columns  (over  or,  exten- 
s i cm  ’ield  containing  ail  of  the  above  c - ' s ) i .e. , r^c^  + -r  .c^  + 
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Where  c4,  c,, c are  columns  of  £,  and  a.,  a,,..,,  are  seal- 

S J »l  1 J 

ers  which  may  be  elements  of  the  extension  field.  Thus  it'  c is  the 

m 

dependent  column,  delete  c from  Ui  and  R from  p.  Add  <*R  to  row  R . , 

mTm  1 m 1 

ajRjjj  to  row  R^,  etc..  Repeat  this  process  until  all  linear  dependen- 
cies among  the  columns  of  ^ have  been  deleted.  Denote  the  resulting 
equations  by  £ p = <5  (3.21) 

Here  6 is  a vector  of  the  extension  field  elements,  £ has  a left 
inverse  given  by,  £ "L  = (£t  £)”'  £t,  (3.22) 

Hence,  one  can  solve  for  p as,  p = £ _L  6,  (3.23) 

which  is  a vector  of  elements  in  the  extension  field.  Finally,  one 
desires  to  compute  the  g^'s  from  p.  One  can  express  p as  £ = eP  where 
6 is  a known  matrix  of  extension  field  elements.  There  is  also  an  add- 
itional constraint  that  the  g^s  lie  In  GF(p),  even  though  b and  p are 
composed  of  extension  field  elements.  Since  every  extension  field  ele- 
ment can  be  uniquely  represented  as  a mth  order  polynomial  in  ai  indeter- 
minant a with  coefficients  in  GF(p),  ^ Equation  (3.23)  can  be  written 

m m 

I [e  rV  = l [q  j«r 
Y=0  y r=0 

where  and  qr  are  the  coefficients  of  the  polynomial. 

Upon  equating  coefficients  of  like  powers  of  a»  one  obtains  a set 
of  simultaneous  equations  in  GF(p)  that  are  to  be  solved  for  the  9^'s. 

These  equations  can  be  solved  by  any  standard  solution  techniques 
in  GF(p).  In  particular,  in  GF(2),  one  can  set  up  a Boolean  express- 
ion, which  after  simplification,  reduces  to  a list  of  all  possible  sets 
of  component  values  consistent,  with  the  specified  data. 


i 
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In  the  following  chapter,  several  examples  of  LSC  fault  analysis 
are  considered  In  which  the  technique  described  above  Is  used  for  solv- 
ing the  set  of  simultaneous  equations  arising  In  the  fault  analysis. 


Figure  4.1:  A LSC  for  Example  4.1 


Assume  that  the  zero  state  impulse  response  for  the  LSC  in 

figure  4.1  is  measured  and  is  given  by,  hQ  = 1 , hj  = 1 , = 1 

It  is  desired  to  compute  all  possible  values  of  the  gains  g1 , g?, 
g3>  g^  which  are  compatible  with  the  zero  state  impulse  response  data 
given  above. 

Solution:  The  connection  equations  for  this  circuit  are  as  follows. 
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0 


5 


I 


Substituting  the  matrices  L...  G and  F(D)  obtained  from  (4.1)  and  (4.2) 

* J 

into  (3.6)  and  then  using  the  term  expansion  algorithm,  one  obtains. 


r 


Equation  (4.3)  Is  of  the  form 
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t p = <5 

where  zero  entries  In  y and  in  the  corresponding  rows  of  p are  deleted. 
The  solution  of  the  set  of  Equation  (4.3)  is  carried  out  in  (GF(2  )). 
Choose  the  elements  e^o,  e2=a  , e3=a  , e4~a  , e^a  , eg=a  from 
GF(23)  where  a Is  an  indetermlnant.  Use  (2.25)  to  compute  S(e.|)  for 
1 * 1 , 2,  3,  4,  5,  6 from  the  zero  state  impulse  response  data.  They 
are  given  as  in  (4.4) 

S(e| ) = aS 
S(e2)  = a3 
S(e3)  = a2 

s(e4^  = 

3(e5)  - a 

S(e6)  * a4  (4.4) 


The  set  of  simultaneous  equations  (4.3)  are. 


1 

1 

S(e] ) 

s(e]) 

S^) 

i s(ei) 

9] 

92 

- — 
Sfe,) 

ei 

ei 

el 

ei  4 

9293 

9294 

— 

9193 

• 

919294 

1 

1 

S('6> 

S(e6) 

S(e,) 

' s<'6> 

91929394 

S(e6) 

e6 

e6 

e6 

e6  e6 

— 

— < — 

— _ 

Substituting  (4.4)  into  (4.5) 
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In  (4.6)  the  matrix,  has  only  three  linearly  independent  columns 
and  can  be  reduced  to  4 by  using  the  equalities. 


(4.-0 


» 


The  equation  (4,7)  is  of  the  form. 
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where. 


or. 


1 1 

0 C 

C 0 


0 

1 

0 


Ijj  D K 6 


+ 9, 


9293  + 9294  + 9193  + 91929394 
919294  + 91929394 


0 

1 

0 


0 

1 

0 


0 

0 

i 


g 

9-> 


9293 

g294 


o.q„ 

' i * j 


•1 


Now  multiplying  both  sides  of  (4.7)  by  one  obtains. 


919294 

[^l929394j 


1 

“ ‘ 

91  + 92 

1 

-L  , 

9.93  ♦ 9294  ♦ g,g3  + g192s3g4 

K ii  iS  r 

1 

Wa  + 91929394 

0 

(4.8) 


Note  that  «il  coefficients  are  in  GF(2).  so  the  g^'s  can  be  solved  for 
directly  by  expanding  the  Boolean  equation,  where  in  GF(2)  is  inter- 
preted as  an  "EXCLUSIVE  OR"  operation  in  Boolean  Algebra. 


(9,+92)  (9293+9294+9193+919293s4)  (919294+91929394)  = 1 


[t  .0) 


Assume  that  the  measured  zero  state  impulse  response  Is 
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{ h } ={  1, 1,0,1, 0,0, M } 


Solution  : matrices  are  same  as  given  in  (4.1)  and 


[ (G  (D)  3 o 0 


0 1/0"  0 


(4.12) 


Substituting  the  matrices  L^. , G and  F(D)  obtained  from  (4.1),  (4.12) 
into  (3.6)  and  then  using  the  term  expansion  algorithm  , one  obtains 


the  set  of  equations. 


1 1 

D 


M 

o 


[S'D)] 


(4.13)  i..  jf  the  form. 


p = 6 


i ^4 


W4 

w3;j4 


(4.13) 


where  zero  entries  in  'i'  anc  in  corresponding  rows  of  p are  dtleted. 

3 

Attempting  to  use  the  e’s  from  GK(?  ),  one  obtains. 


f. 
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S(e,)  = S(a)  = a' 


SCe2)  = s(a2)  = a6 


S(e3)  = S(a3)  = a5/0  i.e.,  undefined 


S(e4)  = S(a4)  = a5 


S(e5)  = S(a5)  = a3/0 


i.e.,  undefined 


S(eg)  = S(a6)  = a6/0  i.e.,  undefined 


Note  that  e3  = a3,  e5  = a5,  efi  = a6  are  poles  of  S(e)  and  hence 

these  three  elements  can  not  be  used  to  set  up  the  set  of  simultaneous 

equations  (2.26).  As  such  one  can  write  three  fault  analysis  equations 

in  GF(23)  which  are  not  sufficient  to  solve  for  the  5^'s.  However,  S(e^) 

is  defined  in  GF(2  j.  Tables  of  multiplication  and  addition  for  GF (2  ) 

15 


are  given  in  reference 


Equation  (2.25)  can  now  be  used  to  compute  $(e^}  from  the  zero 


state  impulse  response  data.  These  are  given  by, 

11 


S(e-|)  = 5(a) 


S(e2)  = $ (a* 


S(e3)  = S(-.t  ) 


S(e4)  - S(<4)  = a14 


v fl 

e5' 


5 ( c 


= 1 


S(efi)  = S(a6)  = a12 


Ke,)  = S(a') 


(4.14) 


Upon  expanding  (4.13)  for  each  ei  from  GF(24)  and  substituting  in  (4.14) 
and  in  (4.13),  the  set  of  simultaneous  equations  obtained  are. 


r 


1 

13 

8 

8 

10 

12 

a 

a 

a 

a 

a 

a 

1 

n 

5 

9 

a 

a 

a 

a 

a 

c. 

1 

9 

12 

12 

3 

6 

a 

a 

a 

a 

U 

a 

1 

7 

2 

2 

10 

3 

a 

a 

a 

a 

a 

a 

1 

5 

a 

i 

i 

10 

a 

i 

a 

1 

3 

9 

9 

6 

12 

a 

a 

a 

a 

a 

a 

1 

2 

2 

9 

a 

a 

a 

a 

a 

a 

*hg3 

gl9294 


^l92g394j  £ J 


Equation  (4.15)  is  of  the  form, 

^ p - 6 

After  deleting  linearly  dependent  columns  of  y and  then  using  the  tech- 
nique described  in  Example  4.1,  one  obtains, 

f~9~  8 [I  9 

a Q-.  ■*  a n-a.a-0 . a 

*1 

- = g-|  + ag2  1 + a 

7 2 ~ 

u 9-|92^4  + n ^1929394  9 

a10g1  + a11g]93  + ;i791929394  <>14  (4.16 


(4.16) 


Equation  (4.16)  yields  a unique  solution, 

91  = 1 ’ g2  = 1 ’ 93  £ 1 » = 0 • (4.17) 

One  can  check  this  solution  by  finding  an  impulse  response  of  the 
LSC  in  Figure  4.2  by  using  the  values  of  gain  in  (4.15). 

The  important  point  here  is  that  with  the  additional  delays,  the 
internal  component  paramfters  can  be  determined  exactly  from  the  imt  1 \p 
response  whereas  without  the  delay  the  best  that  one  can  do  is  to  obtain 
a list  of  three  possible  'aults. 
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CHAPTER  V 

FAULT  ANALYSIS  - AFFINE  SEQUENTIAL  CIRCUITS 
5.1  Introduction: 

Mathematically,  a "two-sided"  Affine  Sequential  Circuit  (ASC)  over 
a finite  field  is  represented  by  a set  of  difference  equations. 


Vi  ■ flxk  * Bv 

yk  = Cxk  * Juk  * wk. 


< k < 


(5.1) 


xo  * 


where  {w^}  is  a sequence  of  constants  ever  a finite  field  and  w^  is 
some  constant  at  time  k.  Other  terms  in  (5.1)  are  as  defined  in  Chap- 
ter II. 

In  cases  where  ASC's  are  not  given  in  the  form  shown  in  (5.1),  a 
change  of  variables  and  some  manipulations  will  yield  the  desired  form 
(5.1).  Hence,  without  loss  of  generality,  (5.1)  can  be  considered  as  a 
standard  form  for  the  purpose  of  fault  analysis  in  Affine  Sequential 
Circuits. 

For  Affine  Sequential  Circuits  the  present  value  of  the  output 
depends  riot  only  on  the  present  input  value  and  the  present  state  (as  is 
the  case  with  LSC's)  but  also  cn  some  constants  defined  over  a finite 
field.  Since  the  mathematical  representation  of  ASC's  and  LSC's  only 
differ  by  some  constant  value  in  the  output  equation,  the  existence  and 
uniqueness  theory  for  the  solution  of  (5.1)  is  similar  to  that  discussed 
in  Chapter  II  for  ISC's  and  will  not  be  repeated  here. 

5.2  Component  Connection  Model : 

Since  Affine  Sequential  Circuits  have  affine  components,  each 
component  of  an  ASC  is  characterized  by  the  D-transfer  equation. 
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bi(D)  = [(gi)fi(D)J  a.(D)  «■  fq.(D)],  i = 1,2 n (5.2) 

where,  a^D)  is  the  D-transform  of  the  input  sequence  to  the  i^ 
component. 

b.j{D)  is  the  D-transform  of  the  ouput  sequence  to  the  itfl 
component, 

f^'Q)  represents  the  dynamics  of  the  ith  component, 
g^,  a scalar  of  the  underlying  finite  field,  is  the  gain  of  the 
linear  part  of  i^  component. 

Let  n • be  a bias  sequence  added  to  the  ith  component.  The  D-trans- 
form  of  this  bias  sequence,  n^(D)  is  given  by, 

ni(D' = ( D^r  )ni 

qH(D)  is,  then,  a product  of  r^(D)  and  whatever  part  of  the  com- 
ponent dynamics  the  bias  signal  passes  through. 

For  the  purpose  cf  fault  analysis,  it  is  assumed  that  the  component 
dynamics  f - (D)  remain  constant  and  all  faults  manifest  themselves  as 
changes  in  g^  and  r,^.  Such  a fault  model  includes  " stuck-or,-one"  faults 
which  are  not  included  in  the  case  of  LSC's. 

One  may  note  that  an  /'SC  over  GF(?)  can  be  viewed  in  either  of  two 
ways;  first  as  an  LSC  in  which  a constant  bias  source  is  introduced  and 
secondly  as  an  LSC  into  which  NOT  gate  has  been  inserted.  These  two 
viewpoints  are  equivalent  since  one  can  construct  a NOT  gate  with  a 
bias  adder  and  conversely. 

For  rotational  simplicity,  the  n scalar  eguations,  (5.4),  can  *• 
combined  into  a single  matrix  equation, 

b(0)  --  [ (G)F(n)]  a(P)  - [Q(n)l  (5.5) 


where 


b(D)  co1umn(b^  (D) ) , 
e (D)  = tolumn(a . (D) ) , 
G = dia9.(gi )» 

F(D)  = diag.(f -(D)), 
Q(D)  = diag . (q . (D) ) . 


Example  5.1: 

This  example  illustrates  the  formulation  of  Equation  (5.5)  for  the 
component  shown  in  Figure  5.1. 


Figure  5.1:  Figure  for  Example  5.1 


Let  G = 1 


Then 


bf  P' 


:m  ♦ < &><  i)  , 


As  in  the  previous  discussion  of  the  component  connection  model, 
the  connection  structure  ot  an  Affine  Sequential  Circuit  is  described 
by  the  algebraic  constraints. 


a(D) 

~ n 

L12 

b(D) 

y(o) 

L?1 

L2? 

u(D) 

(5.6) 


t 


Simultaneous  solution  of  (5,t0  and  (5.6)  yields. 
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y<D)  = [ l22  + L2,(I  - (G)F(O)  Ln  f1  (G)F(D)L12  ] u(D) 

+ [ L21(I  - (G)r(D)  L^)'1  (Q(D) ) ] (5.7) 

Denote  S(P)  = [L22  + S_21  (I  - (G)F(D)Ln  )‘]  (G)F(D)L1?]  and 

T(D)  « [L21(l  - (G)F(D)L11)‘'I(0(D))]  (5.8) 

Then  (5.7)  becomes, 

y(D)  * S(D)  u(0)  + T(D)  (5.9) 

From  (5.8)  it  is  seen  that  the  D-transform  of  a transfer  function, 
S(D),  is  the  linear  part  of  the  ASC. 

5.3  Fault  Analysis: 

Fault  analysis  in  Affine  Sequential  Circuits  involves  developing 
techniques  for  determining  G and  n.  The  unknown,  G,  can  be  solved  for 
by  considering  the  linear  part  of  an  Affine  Sequential  Circuit. 

If  one  can  determine  the  impulse  response  of  the  linear  part  of  an 
ASC,  Equation  (3.15)  may  be  set  up  for  several  e's  in  the  algebraic 
closure  of  a finite  field  and  then  can  solve  for  the  unknown  g^'s  using 
the  algorithm  described  in  Chapter  III.  The  impulse  response  of  the 
linear  part  of  an  ASC  is  determined  as  follows. 

Observe  the  response  of  'he  ASC  due  to  an  input  sequence  lu^l  = {0}. 
When  this  response  becomes  periodic,  apply  an  impulse  sequence  to  the 
ASC  under  study.  The  impulse  response  of  the  linear  part  of  an  ASC  is 
then  qiven  by  the  difference  of  two  responses  i.e.,  response  obtained 
after  and  before  the  impulse  sequence  is  applied. 

After  the  g.'s  are  determined,  it  remains  to  determine  the  i^'s 
Taking  the  D-transform  of  Equation  (5.1)  and  then  comparing  with 


(5.9) , one  obtains. 


55 

W(D)  = T(D)  . (5.10) 

where  W(D)  is  the  !)-transform  of  the  sequence  {wkK 

If  can  be  determined  from  the  measurement  on  an  ASC,  n can  be 
determined  from  (5.10).  The  sequence  {wk'  is  determined  from  the 
measurements  on  an  ASC  as  follows: 

Let  an  input  sequence  {u^}  be  an  impulse  sequence,  then  Is  follows 
from  (5.1)  that 

h'k  = CAkx  + CAk-1B  + wk,  k > 0 (5.10A) 

where  hk  is  the  impulse  response  of  an  Affine  Sequential  Circuit.  Again 
using  the  controllability  criteria  as  discussed  in  Chapter  III,  one  can 
write,  k 

CAkx  = l CAkR  b (S.ll ) 

“ jSi  -J 

where  b ^ , kQ  < -j  < 0 is  a sequence  of  inputs  that  drives  the  minimal 
realization  of  an  ASC  from  the  zero  state  at  time  k = ko  to  x_  at  time 
k = 0.  Such  a sequence  exists  for  a minimal  realization  if  kQ  < ~6, 
where  f is  the  dimension  of  the  minimal  realization.  One  has  the 
equality, 

[CAk'1B]  = hk  (5.12) 

where  hk  is  the  impulse  response  (at  time  k)  of  the  linear  part  of  an 
ASC. 

Substituting  (5.11)  and  (5.12)  into  (5.10)  yields, 
ko  k . 

h'  = \ CA b . + h,  + w. , k > 0.  (5.13) 

A j_  "j  " J A * 

From  (5,13)  one  can  solve  for  w,  (and  b ).  No  attempt  is  made 

here  to  develop  any  technique  for  solving  wk  from  (5.13).  Instead  an 


I 


t 


I 
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assumption  is  made  that  xQ  = x - 0 which  simplifies  the  expression 
(5.10),  yielding 

hk  = hk  * wk 
or 

»k  * h'  - hk  (5.)4) 

With  the  above  theory,  one  can  formulate  the  fault  analysis  algo- 
rithm on  ASC's.  The  procedure  consists  of  the  following  steps: 

(i)  Solve  for  gp  from  the  impulse  response  of  the  linear  part  of  an 
ASC  and  then  use  steps  (i)  through  step  (iv)  given  in  Chapter  III  for 
the  purpose  of  fault  analysis  on  LSC's; 

(ii)  Substitute  the  value  of  g„.  's  (obtained  from  step  (ij)  in  T(D); 

(iii)  Measure  the  impulse  response  of  the  Affine  Sequential  Circuit  and 
using  (5.14)  obtain  (wk>; 

(iv)  Obtain  the  D-trar.sform  of  the  sequence  {wfc}  i.e.v  obtain  W(D); 

(v)  Using  equality  (5.10),  solve  for  o,  since  n is  contained  in  T(D). 

5.4  Examples: 

In  this  section  two  examples  are  presented.  Both  examples  illus- 
trate the  fault  analysis  algorithm  discussed  in  the  section  5.3.  The 
second  example  also  illustrates  how  stuck-on-one  faults  are  modeled 
using  Affine  Sequential  Circuits. 

Example  5.2: 


Consider  the  ASC  shown  in  Figure  5.2  defined  over  GF(2). 


1 

i 


Figure  5.2  ; Figure  for  Example  5.2 
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Assume  that  the  Impulse  response  of  the  linear  part  of  an  ASC  Is 

measured  and  Is  given  by  hQ  ■ 0,  h1  « 1 , hg  * 1 > h3  * 1 , » The 

impulse  response  of  the  Affine  Sequential  Circuit  is  measured  and  is 

given  by  h^  * 0,  hj  « 0,  !»£  * 1.  h^  « 0,  h^  *»  1,  « 0,.. It  is 

desired  to  compute  the  values  of  and  that  are  compatible  with  the 
given  data. 

Solution: 

The  connection  equations  for  this  circuit  are  as  follows: 


and 

b,<0)  = 19,/Dj  a,(D)  ♦ (j^)  [1/D]  ni  (5.16) 

5(D)  • [0  ♦ 1(1  - g-j/D  • l)"1  g}/D  • 1] 

- (1  - 9/D)"1  g^D  (5.17) 

T(D)  1 (g?T>  HO-  g-j/0  ’ D’^l/D)] 

« 0 - 91/0)“1(giT)  (5.18) 

Calculating  g, : 

« 

2 

Choose  an  element  e1  * a from  GF(2  ) where  a is  an  indetermlnant. 
Then  use  (2.25)  to  compute  S(e^)  - S{o)  from  the  given  impulse  response 
of  the  linear  part  of  an  ASC. 

From  (2  25)  one  obtains, 

SU)*'  = + 1 5.19) 

S(a)o2  = b , + a + 1 


(5.20) 


Simultaneous  solution  of  (5.19)  and  (5.20)  yields 
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S(a)  " a . 


(5.21) 


Since  S(e.)  = S(D)  n = S(a),  one  can  replace  D by  a in  (5.17) 

I u~e.j  =a 


and  obtain. 


5(a)  - (1  - ^ 


Simultaneous  solution  of  (5.21)  and  (5.22)  yields 

g1  = 1. 

Substitute  1 1 in  (5.18)  to  obtain, 

T(C)  = (1  - ^)-1(oTt) 

* {b~r]  Td^TT 


(5.22) 


(1+0)' 


T i P.C?  — *1  = in  flP  f ? ll 


(5.23) 


1+0 


Obtaining  {wfe}: 

from  (5.14)  one  obtains, 

" 0 - 0 r 0 

0 0 0 

w,  - h,  - h;  - 1 - 0 = 1 

*2  ‘ h2  ' h2  = 1 ‘ 1 = 0 

*3  “ h3  ‘ h3  ' 1 - 0 = 1 

so,  i wk  1 = {0,  1,  0,  1, 1 , such  that  «o  = 0. 

W(D)  - D{wk»  - 'G-n°  + 1-D'1  + 0-D"2  + 1-D~2  + 

w(0)  = - K 

1+D^ 


(5.24) 


0 


solving  for  h: 

Substituting  (5.23)  and  (5.24)  in  (5.10),  i.e., 

WfO)  = 1(0)^ 

D . 0 . 

I.e.,  — -y  - y n, 

1+u  1+D"  1 


whicn  yields  - 1.  (5.25 

It  is  interesting  to  note  fro«n  Figure  5.2  that  is  exactly  oppo- 
site to  a^  which  is  a characteristic  of  the  NOT  gate.  So  ASC’s  also 
permit  modeling  of  NOT  gate  for  the  purpose  of  fault  analysis. 

Example  5.3: 

Consider  the  A$C  shown  in  Figure  5.2.  Assume  that  the  impulse 
•espcnse  of  the  linear  part  of  the  ASC  is  measured  and  is  given  by 

= 0,  hj  * 0,  = 0,  The  impulse  response  of  the  Affine 

■^cULiential  C-.rcuit  is  measured  and  is  given  by  - 0,  h.j  ' 1,  lij  ••  1, 
r:  - 1,  It  is  desired  to  compute  the  values  of  g.  and 

3 I 1 

that  are  compatible  with  the  given  data. 

The  connection  equations  are  the  same  as  those  of  the  previous 
example.  Thai  is,  they  are  given  by  (5.15),  (5.16),  (5.17)  and  (5.18). 
Calculating  g, : 

9 

Choose  an  element  e,  - 1 frem  GF(2  ).  Use  (2.'25)  and  the  impulse 
response  of  the  linear  part  of  an  ASC  to  compute  S(n).  This  is  oiven 
ty  S{  «)  - (5.2F 

S i nru  1 taneous  solution  of  (5.26)  and  (5  22)  yields  g^  " 0. 


Suust.tute  1 G in  (5.18)  to  obtain, 


(5.27) 


Obtaining 
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so. 


w ■ h - h"  * 0 
0 0 0 


*1  = hl  ’ h1  a 1 
»2  5 ^ - h£  = 1 


{w^}  = (0,  1 , 1 , 1 , 1 , . . , . } such  that  « 0, 


W(D)  = 0{wk)  = {O’D0  + I'D-1  + 1-D‘2  + 1 *D~3 


(5.28) 


Substituting  (5.27)  and  (5.28)  in  (5.10)  yields  * 1 . 

The  above  example  is  an  illustration  of  a "stuck-on-one"  fault. 


CHAPTER  VI 


CONCLUSION 


In  this  work,  a spectral  theory  for  Linear  Sequential  Circuits  has 
been  formulated  and  the  component  connection  theory  has  been  applied  to 
ISC's.  From  the  spectral  theory  and  the  connection  function  of  the  com- 
ponent connection  theory,  a fault  analysis  procedure  for  a Linear  Sequeri 
tial  Circuit  has  been  developed.  This  procedure  parallels  the  multifre- 
quency testing  technique  for  fault  analysis  in  analog  circuits.  It  has 
beer  shown  that  it  is  often  easier  to  do  fault  analysis  in  sequential 


circuits  than  in  combinatorial  circuits. 

This  fault  analysis  algorithm  has  been  extended  to  include  Affine 
Sequential  Circuits.  The  fault  analysis  procedure  for  Affine  Sequential 
Circuits  has  turned  out  to  be  no  harder  than  fault  analysis  in  LSC's 
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example  that  Affine  Sequential  Circuits  can  be  used  for  modeling  nonlin- 
ear faults  such  as  "stuck-on-one"  faults.  Affine  Sequential  Circuits 
should  cover  a broad  range  of  digital  circuits  hence  an  interesting  area 
for  further  research. 

Not  much  work  is  done  in  the  area  of  fault  analysis  in  nonlinear 
analog  or  digital  circuits.  Due  to  the  finiteness  of  states  of  sequential 
circuits,  one  may  possibly  generate  fault  analysis  algorithms  for  nonlin- 
ear sequential  circuits.  Such  algorithms  can  then  be  generalized  to  in- 
clude large  scale  digital  systems. 
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APPENDIX  A 
GALOIS  FIELD  THEORY 

Finite  Fields  : A finite  field  is  a field  with  a finite  num- 
ber cf  elemerts.  These  fields  are  known  as  Galois  Fields  In  honor  of 
French  mathematician  who  first  Investigated  their  properties. 

Let  I denote  the  integral  domain  of  all  integers  and  p be  a prime 
number.  Consider  the  system  I/(p),  whose  elements  are, 

0 « (0  + Bm  ) 
m 

o *nm  i 

FT  = (p-l  + Bm  ) 

where  Bm  = { b|b=kp,  pel,  k=0,l,2,3, } 

If  a and  c are  elements  of  I/(p),  addition  and  multiplication 
are  given  by, 

a + c = a + c 
a c = a c 

It  is  known  that  1/ ( p ) forms  a finite  field.  This  field  is  denoted 
by  GF(p).  For  example,  GF(2)  has  two  elements  They  are  0 and  T.  It's 
addition  and  multiplication  tables  are  given  below  : 

+ 0 T 0 T 

0 o T . ooo 

T o T o T 
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1 
\ ] 


T 


L'O 

Extension  F i e 1 cis  : If  £ is  a fieM  and  F is  a subfield  of  E, 

£ is  an  extension  of  F.  The  relation  of  being  an  extension  of  F is 
denoted  by  F £. 

The  following  theorems  characterize  the  properties  of  an  extension 

1 3 

field.  Their  proofs  are  given  iri  >efernces 

Theorem  1 : Each  Galois  field  contains  a unique  sut  field  with  a 
p-  iire  number  of  elements. 

Theorem  2 : For  every  positive  integer  in  and  a prime  r.u:  her  p, 
there  exists  an  irreducible  polynomial  in  GF(p)  of  degree  m and  these 
fields,  denoted  by  GF(pm),  have  pm  elements.  GF(pm)  is  an  finite 
extension  of  F. 

Kronecker  Theorem  : If  f(x)  is  a polynomial  with  coc-Tf icier.ts  in 
a field  F,  there  exists  an  tension  F of  F si«  h '.’'at  for  - ' -,e  x in  F, 
f(xQ)  - 0. 

The  algorithm  for  generating  elements  of  GF(pm)  is  described  below: 

Pick  an  irreducible  polynomial  P(x)  of  degree  m.  Introduce  a ■.  ?w 


symbol.-!,  and  as'-ome  tJ‘  :)  - 0.  Then  0,1,. 
pnl  elements  of  Gf  ( pnt ) su.h  that 


p -p  will  be  <;  set  of 


i) 


r 

v 


i i ) -i 1 u'1 


■ , i .ir  . 
v ■ ) '3d  p -p+1 


APPENDIX  B 


DEFINITION  OF  THE  D-TRANSFQRK 

Let  I denote  natural  integers  and  G be  the  set  of  all  sequences. 
Let  g : I -*•  GF(pm)  for  all  g c G. 

The  H-transform  of  g(I)  Is  deftned  in  terms  of  formal  power  series 
in  indeteminant  D denoted  by  D T g f 1 } > - 

D{g(1)}  * VgdJO*1 

\ C.flD 

= +...+D3g(-3)  + D2g(-2)  + Dg(-l)  + g(0)  + D"1g(l)+ 

D"2g{2)  .+D"kg(k)  


L 


p j 
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